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1 Introduction 

Generalized Parton Distributions (GPDs) [1] parametrize the non-perturbative 
hadron structure in hard exclusive processes [2]. The measurement of GPDs 
would represent a unique way to access crucial features of hadron structure 
[3, 4]. Relevant experimental efforts to measure GPDs, by exclusive electron 
Deep Inelastic Scattering (DIS), will take place soon and it is urgent to produce 
predictions for these quantities. Recently, calculations have been performed in 
Constituent Quark Models (CQM) [5, 6]. The CQM has a long story of success- 
ful predictions in low-energy studies of the structure of the nucleon. At high 
energies, in order to compare predictions with data, one has to evolve, according 
to perturbative QCD, the leading twist component of the structure functions 
obtained at the low- momentum scale, the "hadronic scale" /ig, associated with 
the model. Such a procedure, already addressed in [7], has proven successful in 
describing the gross features of standard Parton Distribution Functions (PDFs) 
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(see, e.g., [8]). Similar expectations motivated the study of GPDs in Ref. [5], 
where a simple formalism has been proposed to calculate the quark contribution 
to GPDs. The procedure of Rcf. [5] has been extended and generalized in Ref. 
[9]. As a matter of fact, the approach of Ref. [5], when applied in the forward 
case, has been proven to be able to reproduce the gross features of PDFs [8] 
but, in order to achieve a better agreement with data, it has to be improved. 
In a series of papers, it has been shown that DIS data are consistent with a low 
energy scenario dominated by composite constituent quarks of the nucleon [10], 
defined trhough a scheme suggested by Altarclli, Cabibbo, Maiani and Petronzio 
(ACMP) [11], properly updated. The same idea has been recently applied to 
show the evidence of complex objects inside the proton [12], analyzing data of 
electron scattering off the proton. In this talk, we review the main findings of 
Ref. [9], where the same idea has been applied to obtain realistic predictions 
for GPDs and, at the same time, to explore kinematical regions, not accessible 
with the approach of Ref. [5]. In particular, the evaluation of the sea quark 
contribution becomes possible, so that GPDs can be calculated in their full 
range of definition. Such an achievement will permit to estimate the relevant 
cross-sections, providing us with a tool for planning future experiments. 

2 GPDs in a constituent quark scenario 

Wc arc interested in hard exclusive processes, such as Deeply Virtual Compton 
Scattering (DVCS). We use the formahsm of Ref. [4]. Let us think now to 
a nucleon target, with initial and final momenta P and P', respectively, with 
A = P' — P being the momentum transfer. The main quantity we discuss 
is the GPD H{x,^,A'^). In a system of coordinates where the virtual photon 
4-momentum, = {qo,q), and P = (P + P')^ /2 are collinear along z, the 
^ variable, the so called "skewedness" , parameterizing the asymmetry of the 
process, is ^ = — n • A, where n is a light-like 4-vector with n ■ P = 1. With 
this choice, GPDs describe the amplitude for finding a quark with momentum 
fraction x + (in the Infinite Momentum Frame) in a nucleon with mo- 
mentum (1 -|- ^/2)P and replacing it back into the nucleon with a momentum 
transfer A**. Besides, when the quark longitudinal momentum fraction x of the 
average nucleon momentum P is less than — (the so-called negative DGLAP 
region), GPDs describe antiquarks; when it is larger than ^/2 (the positive 
DGLAP region), they describe quarks; when it is between — ^/2 and ^/2 (the 
ERBL region), they describe qq pairs. One should keep in mind that, besides 
the variables x,£, and A^, GPDs depend, as the standard PDFs, on the mo- 
mentum scale at which they are measured. For an easy presentation, such 
a dependence will be omitted. In [5], the Impulse Approximation (lA) expres- 
sion for Hq{x, ^, A^), suitable to perform CQM calculations, has been obtained. 
In particular, it has been shown that, taking into account the quarks degrees 
of freedom only, considering a process with A^ <^ M^, in a non relativistic 
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Figure 1: The DVCS process in the present approach. 



framework one has 

J dk s(^x + ^-^^ ng{k,k + A) 

j dk 5(x + ^-^^ J e'«^^+^)^^-^^^")p,(f,r')dfdr' ,(1) 

where the one-body non diagonal charge density p{f, f) and the one-body non- 
diagonal momentum distribution hq{k, k + A) have been introduced. The above 
equation allows the calculation of Hq{x,^, A"^) in any CQM, and it naturally 
verifies the two crucial constraints of the GPDs [5]. With respect to Eq. (1), a 
few caveats are necessary. 

i) If use is made of a CQM, containing only constituent quarks only the 
DGLAP region can be explored. In order to introduce the study of the ERBL 
region, so that observables can be calculated, the model has to be enriched. 

ii) in Eq. (1), the x variable for the valence quarks is not defined in its natural 
support. Several prescriptions have been proposed in the past to overcome such 
a difficulty in the standard PDFs case [7, 8]. 

iii) The Constituent Quarks arc assumed to be point-like. 

The two issues i) and ii) will be now discussed, by relaxing the condition iii) 
and allowing for a composite structure of the constituent quark 

The procedure described in the previous section, when applied in the stan- 
dard forward case, has been proven to be able to reproduce the gross features 
of PDFs [8]. In order to acliieve a better agreement with data, the approach 
has to be improved. In a series of previous papers, it has been shown that DIS 
data are consistent with a low-energy simple picture of the constituent quark as 
a complex system of point-like partons. [10], updating the ACMP scenario [11]. 

Following the same idea, we describe here a model for the reaction mecha- 
nism of an off-forward process, such as DVCS. As a result, a convolution formula 
giving the proton Hg GPD in terms of a constituent quark off-forward momen- 
tum distribution, Hq^, and of a GPD of the constituent quark go itself, Hg^q, 
will be derived. It is assumed that the hard scattering with the virtual photon 
takes place on a parton of a spin 1/2 target, made of complex constituents. The 
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scenario we are thinking to is depicted in Fig. 1 for the special case of DVCS. 
In addition to the kinematical variables, x and ^, already defined, one needs few 
more ones to describe the process. In particular, x' and for the "internal" 
target, i.e., the constituent quark, have to be introduced [9]. Using lA, a stan- 
dard procedure can applied and a convolution formula, is readily obtained (see 
Ref. [9] for details): 

go ^ ^ ^ 

where Hq^{z,^,^) is to be evaluated in a given CQM, according to Eq. (1), 
for go = Wo or do, while |, A^) is the constituent quark GPD, which 

is still to be modelled. We can start modelling this quantity thinking first of 
all to its forward limit, where "constituent quark parton distributions" have 
to be recovered. As wc said, in a scries of papers [10] a simple picture of the 
constituent quark as a complex system of point-like partons has been proposed, 
re-taking the ACMP scenario, [11]. 

Let us recall the main features of that idea. The constituent quarks are 
complex objects whose structure functions are described by a set of functions 
4>qgq{x) that specify the number of point-like partons of type q which are present 
in the constituent of type go, with fraction x of its total momentum. Wc will 
call these functions the structure functions of the constituent quark. The func- 
tions describing the nucleon parton distributions are expressed in terms of the 
independent (f)qQq{x) and of the constituent density distributions (go = uo,do) 
as, 

q{x, Q') =Y,f -Qo{z, Q')^q,q (-, Q^) , (3) 

where q labels the various partons, i.e., valence quarks {uv,dv), sea quarks 
{us,ds,s), sea antiquarks {u,d. s) and giuons g. The different types and func- 
tional forms of the structure functions of the constituent quarks are derived 
from three very natural assumptions [11]: i) The point-like partons are QCD 
degrees o freedom, i.e. quarks, antiquarks and giuons; ii) Regge behavior for 
X ^ and duality ideas; iii) invariance under charge conjugation and isospin. 
The last assumption of the approach relates to the choice of the scale at which 
the constituent quark structure is defined. We choose /i§ — 0.34 GeV^ [8, 14]. 
This hypothesis fixes all but one the parameters of the approach. The remain- 
ing paramter is fixed according to the value of F2 at low x[ll], and its value is 
chosen again according to [14]. The values of the parameters obtained are listed 
in [10]. The other ingredients appearing in Eq. (3), i.e., the density distribu- 
tions for each constituent quark, are defined according to Eq. (1). Now we have 
to generalize this scenario to describe off-forward phenomena. Of course, the 
forward limit of our GPDs formula, Eq. (2), has to be given by Eq. (3): this 
gives if,„,(f,|,A2) [9]. 
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Now the off-forward behaviour of the Constituent Quark GPDs has to be 
modelled. This can be done in a natural way by using the "a-Double Distribu- 
tions" (DD's) language proposed by Radyushkin [13]. DD's, $(i,a, A^), rep- 
resent an alternative parametrization of the matrix elements describing DVCS 
and hard exclusive clcctroproduction processes, with respect to the one based 
on GPDs. The DD's do not depend on the skewedness parameter ^; rather, they 
describe how the total, P, and transfer, A, momenta are shared between the 
interacting and final partons, by means of the variables x and a, respectively. 
Hqgq for the constituent quark target, is related to the a-DD's, which we call 
$5(55(5, a, A^) for the constituent quark, in the following way: 

Hg^g{x,^,A'^) = J ^dx J ^^6(^x + ^a-x^^g^q{x,a,A'^)da. (4) 

In [13], a factorizcd ansatz is suggested for the DD's: 

^g,g{x,a,A^) = hg{£,a,A'')<S?g„giS:)Fg,iA'') (5) 

where the a dependent term, hq(x, a, A'^), has the character of a mesonic am- 
plitude, ^qgq{x) represents the forward density and, eventually, Fqg{A'^) the 
constituent quark f.f. In the following, we will assume for the constituent quark 
GPD the above factorizcd form, so that we need to model the three functions 
appearing in Eq. (5), according to the description of the reaction mechanism 
we have in mind. For the amplitude hq, use will be made of one of the sim- 
ple normalized forms suggested in [13]. Besides, in our approach the forward 
densities $gog(i) have to be given by the standard $ functions of the ACMP 
approach (see Ref. [9] for details). Eventually, as a f.f. we will take a monopole 
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Figure 3: The NS H GPD for valence u-quarks at C = 0.2 and = -0.5 
GeV^, at the momentum scale of the model, /Iq = 0.34 GeV^ (dashed), and 
after NLO-QCD evolution up to = 10 GeV^ (full). 



form corresponding to a constituent quark size rq ~ 0.3/to, supported by the 
analysis of [12]. By using these ingredients, Eq. (4) can be calculated and the 
result, when cast into Eq. (2), yields the GPD in our ACMP approach. 

Results for the GPD Hq of the proton have been obtained calculating Eq. 
(2). Hq^q, has been calculated according to the model discussed above, and Hq^ 
has been evaluated according to Eq. (1), using, as an illustration, the wave 
functions of the Isgur and Karl (IK) model [15]. A prescription introduced in 
[9] has been used to correct the poor-support problem, addressed in Section 2. 

Results are shown in Fig. 2. A relevant tail of the valence quarks contribu- 
tion in the ERBL region is found, in agreement with other estimates [2]. The 
knowledge of GPDs in the ERBL region is necessary for the calculation of cross- 
sections. The ERBL region is accessed here, with respect to the calculation of 
Ref. [5], thanks to the constituent structure which has been introduced. 

The results shown so far have to be related to the hadronic scale ij,q=0.34: 
GcV^. In Fig. 3 the NLO QCD-evolution of the Non Singlet, valence u distri- 
bution, is shown. In evolving the results, the approach of Ref. [16] has been 
applied and a code provided by A. Preund has been used. The evolution clearly 
shows a strong enhancement of the ERBL region. 

3 Conclusions 

The aim of the talk has been to describe composite constituent quarks [11, 10] in 
studies of GPDs [9] . This is the continuation of an effort to construct a scheme 
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which describes hadrons in different kinematical and dynamical scenarios. We 
have developed a formalism which expresses the hadronic GPDs in terms of 
constituent quarks GPDs by means of appropriate convolutions. Looking at 
our results, we discovered that such a scheme transforms a hadronic model, 
in whose original description only valence quarks appear, into one containing 
all kinds of partons. Moreover, the starting model produces no structure in 
the ERBL region, while after the structure of the constituent quark has been 
incorporated, it does. The completeness of the a;-range, for the allowed and 
^, of the present description, permits the calculation of cross-sections in a wide 
kinematical range and it can be used therefore to guide future experiments. 
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